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Abstract 

A differential system [A] : Y' = AY, with A G Mat(n, k) is said to be in reduced form if A G g(fc) where g 
is the Lie algebra of the differential Galois group G of [A] . 

In this article, we give a constructive criterion for a system to be in reduced form. When G is reductive 
and unimodular, the system [A] is in reduced form if and only if all of its invariants (rational solutions of 
appropriate symmetric powers) have constant coefficients (instead of rational functions). When G is non- 
reductive, we give a similar characterization via the semi-invariants of G. In the reductive case, we propose 
a decision procedure for putting the system into reduced form which, in turn, gives a constructive proof of 
the classical Kolchin-Kovacic reduction theorem. 
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1. Introduction 

The direct problem in differential Galois theory is, given a linear differential system, how to compute its 
differential Galois group. Though there exist several theoretical procedures to achieve this task (Hrushovski 
[Hru02], Compoint-Singer [CS99, vdPS03], van der Hoeven [vdH07]), it is still far from being practical. 
If we now ask directly for the Lie algebra, not much is known. Nevertheless one can observe that if Y' = AY 
is a differential system with coefficients in k = C(z) for example, then one can give a "supset" for the Lie 
algebra g of the differential Galois group G. Indeed, Kolchin showed that g is included in any Lie algebra 
f) such that A G i)(k) ([vdPS03], chap. 1). If A G f)(fc), and if we transform this system by a gauge matrix 
P G GL n {k) in such a way that the new system Y' = P[A]Y satisfies P[A] G f)(fc) with () C (), then we 
obtain the finer majoration g C t). The problem is then to find a matrix P such that P[A] G fl(fc). A result 
of Kolchin and Kovacic (proposition 3 page 3) shows that it is always possible to find such P G GL n (k). 
When P[A] G fl(fc), we'll say that the system is in reduced form. But how can one recognize that a system is 
in reduced form? In this paper we give an answer to that question. More precisely, we show (theorem 1) that 
the system is in reduced form if and only if all its semi-invariants can be written in a very particular form: 
they are all the product of an exponential and a constant coefficients vector (instead of general rational coef- 
ficients, i.e. elements of k, as usual). We call them semi-invariants with constant coefficients (see definition 



Email addresses: aparicioOrisc .uni-linz . ac . at (Ainhoa Aparicio-Monfortc*), Elie.CompointOmath.uiiiv-lillel.fr 
(Elie Compoint), Jacques-Arthur.WeilOunilim.fr (Jacques- Arthur Weil) 
* Supported by the Austrian FWF grant Y464-N18 



Preprint submitted to Elsevier 



Sept 2011-May2012 



23). In the case where the Galois group G is reductive and unimodular, this criterion becomes (theorem 2): 
the system Y' = AY is in reduced form if and only if all its invariants have constant coefficients (instead of 
general rational coefficients) . We then propose a procedure which puts a system into reduced form (this also 
gives an algorithmic proof of the Kolchin-Kovacic reduction proposition 3 page 3). 

This work originates in the first author's PhD ([AM10]) where a weaker version was presented. Other ap- 
plications of reduced forms (to intcgrability of hamiltonian systems) are studied in [AMW11, AMW12]. 

The paper is organized as follows. In the first part we recall some basics on differential Galois theory. 
We define, inspired by works of Wei and Norman, the Lie algebra associated to a matrix A £ M n {k) and 
the notion of reduced form. In the second part, we recall standard facts about the tensor constructions 
of a differential module and the representations of the differential Galois group associated with them. In 
the third part we give and prove our main criterion to characterize when a linear differential system is in 
reduced form. The fourth part is devoted to the reductive unimodular case and the algorithmic aspects of 
these criteria. We finish with examples illustrating the results of the paper. 

Acknowledgements. We would like to thank P. Acosta-Humanez, E. Hubert and M.F Singer for fruitful 
discussions on this material. We also thank the referee for subtle and constructive comments. 

2. Reduced Forms of Linear Differential Systems 

2.1. Differential Galois Group 

Let (fc, d) denote a differential field of characteristic zero whose constant field C is assumed to be al- 
gebraically closed. The usual case is k — C(x) with 9 = 4-. We consider a linear differential system 
[A] : Y' = AY with A g Mat (re, k). We refer to the book [vdPS03] for the differential Galois theory used 
here. Let K denote a Picard-Vcssiot extension for [A]. Let G := Aute(K / k) be the differential Galois group 
of [A] and g = Lie(G) its Lie algebra. We also write Lie([A]) for the Lie algebra of the Galois group G of 
Y' = AY. The C-vector space of solutions of [A] in K n is noted V. Remark that the given of a fundamental 
solution matrix U <G GL(n,K) of [A] is equivalent with the choice of a basis (lj)j=i— « of V". 
The differential module associated to [A] is noted Ai = (k n , V^). The connexion Va admits A as its matrix 
in the canonical basis (ej)j = i... n of k n and is defined by V>i(Y) = dY — AY for Y in k n . 

Definition 1. We say that two differential systems Y' = AY and Z' = BZ, with A, B g Mat(n, k) are 
(gauge) equivalent (over k) if there exists a linear change of variables (a gauge transformation) Y = PZ 
with P g GL(n, k) changing [A] to [B], i.e. B = P[A] := p- 1 ■ (A ■ P - d{P)). 

Changing system [A] to an equivalent one is the same as changing basis in the differential module M.. 
Choosing a basis of V yields a faithful representation of G in GL(n, C). Hence there exists a polynomial 
ideal / C CfA^i, . . . , Xij, . . . , X njn , j^-A (where jj^ represents a solution u of u. det(Xij) = I) such that 

G ~ {AI = (rm,j) G GL(n, C) : VP g I , P(m id ) = 0}. 

We call I the ideal of relations of G. Similarly, we obtain a representation of Lie(G) as 

Lie(G) := {N g Mat(n, C) : Id + e N g G(C[e] ) with e ^ and s 2 = 0} 

where Id is the identity matrix, and the C[£r]-points of G are the matrices M £ with coefficients in C[e] which 
satisfy all the equations of the ideal /. 
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2.2. Definition and Existence of Reduced Forms 

Once we have fixed a representation of the differential Galois group G (e.g. by choosing a basis of 
the solution space V of [A]), its Lie algebra g is a C-vector space generated by matrices {Mi, . . . , Md} C 
Mat(n, C). Another choice of basis yields a conjugate representation of g (see remarks below). The set g(k) 
of fc-points of g is g(k) := g <£> k = { /1M1 H h / d Af d , /,£*;}. 

Definition 2. Let A £ Mat(n, fc). Consider the differential system [A] : Y' = AY with Galois group G and 
its Lie algebra g = Xie([y4]). We say that [A] is in reduced form when A £ g(k). 

The following classical result of Kolchin and Kovacic shows that any system admits a reduced form. 

Proposition 3 (Kolchin-Kovacic, [vdPS03] prop 1.31 & cor 1.32). Let k be aC\ field. Consider a differential 
system Y' = AY with A £ Mat(n,fc). Assume that the differential Galois group G (over k) is connected and 
let g be its Lie algebra. Let H C GL(n,C) be a connected algebraic group whose Lie algebra f) satisfies 
A e \){k). Then : 

1. g C f) and G C H. 

2. There exists P £ H{k) such that the equivalent differential system Z' = AZ , with Y = PZ and 
A = P[A], satisfies A e g(k), i.e. A is in reduced form. 

Corollary 4. With the same notations, assume now that G is not connected. There exists P € H{k) such 
that the equivalent differential system Z 1 = AZ , with Y = PZ and A = P[A], satisfies A G g(fc) 

Proof. Let K be a Picard-Vessiot extension of k for Y' = AY. By Galois correspondence (e.g. proposition 
1.34 p.26 in [vdPS03]), if we let k° := fcfl^. we have G ° = Aut d (K/k°). As k° is an algebraic extension 
of k, it is still a Ci-ficld. Pick k° as a base field: then if is a Picard-Vessiot extension of k° whose Galois 
group G° satisfies the conditions of proposition 3. □ 

The assumption that k be a C\ field is used only for the rationality issue (i.e. find a reduced form 
with coefficients in k when G is connected). If we allow a reduced form to have coefficients in an algebraic 
extension of k, then this hypothesis is not needed (we will come back to this and give a short proof of this 
in section 5.3). 

2.3. Wei-Norman Decompositions and Reduced Forms 

Let A = (aij) G Mat(n, k). We consider the C-vector space generated by the n 2 coefficients cii j. Let 
a±,...,a r be a basis of this space; each cijj then admits a unique decomposition a,j = 5Zfe=i TO i^ a fc with 
Tn\j £ C. We thus have, once the basis (ak) is chosen, a unique decomposition A = X)l-=i a kMk where 
M k = (mg) e Mat(», C) 

Definition 5. Let A = (flij) £ Mat(n, k). A Wei-Norman decomposition of A is a decomposition 

r 

A = Y^ a kM k , M k £ Mat(n, C) 

k=l 

where {ai,...,a r } are a basis of Span^^djj)^^! n ). 
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Of course, this decomposition is not unique (it depends on the choice of the basis {ai, . . . , a r }) but its 
dimension r is. Note that the C-vector space generated by the Mj (also of dimension r) depends only on A 

(not on a choice of basis). Indeed, consider an alternative basis {bi, . . . , b r } of Span c ((aij)i J= i r ) such 

that cij = 2i=i Iji^i with Q := (qij)i.j G GL(r, C); this gives an alternative Wei-Norman decomposition 
A = biNi with Ni := Y?j=i QjiMj- It is clear that span c (Mi, . . . , M r ) = span c (A r i, . . . , N r ). 

Definition 6. In the above notations, the Lie algebra associated to the matrix A, denoted Lie(A) is the 
Lie algebra generated by the matrices Mi,...,M r of a Wei-Norman decomposition of A (the vector space 
generated by Mi, . . . , M r and all their iterated Lie brackets). 

Remark 7. The notion of Wei-Norman decomposition borrows its name from Wei and Norman who, in 
[WN63, WN64], use this notion of Lie(A) (also presented there as the Lie algebra generated by all the A(t) 
for t spanning C) to establish their method for solving linear differential systems. 

The above proposition 3 of Kolchin and Kovacic shows that Lie([A]) C Lie(A). The system is in reduced 
form when we reach the equality. 

Definition 8. Y' = AY with A G Mat(n, k). A matrix P G GL(n, k) is called a reduction matrix for [A] if 
P[A] is in reduced form. It is called a partial reduction matrix for [A] if Lie{P[A\) C Lie(A). 

Note that a reduction matrix P may have coefficients in k<~)K. It follows that a reduced form R := P[A] 
may have algebraic coefficients and R G fj(fc). As the change of variables Y = PZ has algebraic coefficients, 
the Galois group may change but not its connected component of the identity, nor its Lie algebra. 

Remark 9. Consider a system Y' = BY in reduced form. Then there exist elements hi G g, a basis {yj} of 
V and a Wei-Norman decomposition B = J2i=i biBi of B such that Bi is the matrix of hi in the basis {y{\. 
Let {iji} be another basis of V , and let Q G GL(n, C) be the matrix of the change of basis from {yi} to {jjj}. 
Then there exists a reduced form Z' = AZ in which the matrices Ai of a Wei-Norman decomposition of A 
appear as matrices in the basis {iji} of elements hi of g. To see this, consider the matrix A = Q~ X BQ. For 
each i, Ai = Q~ x BiQ is the matrix of the element hi of q in the basis {jji}. 

Hence, having a reduced form Y' = AY associated with some a priori given basis {yi} of V (i.e. A = 
Si=i fi ' Ai where the Ai are elements of q represented in the basis {yi}), we may obtain by conjugation 
a reduced form Z' = RZ associated with any choice {iji} of basis of V. When k = C(x), we'll apply this 
result in particular to solutions {iji} from a (local) fundamental solution matrix U determined by the initial 
condition U(zq) = Id at some ordinary point zq G C for A. 

We now prove a simple lemma useful in the proof of our criterion for reduced forms. 

Lemma 10. Let A = ^ i f%A^ be a Wei-Norman decomposition of A. Assume that A admits a constant 
eigenvector w G C n associated with an eigenvalue f G k. Then w is an eigenvector of each Ai. 

Proof. Let w = (ci,... ,c„) ^ 0. Assume, for example, that Cx 7^ 0. Relation Aw = fw gives us fci = 
ctifi with a,j G C, i.e. there exist constants aj such that f = ^2 o-ifi- As Aw = fw, we have J^. fi(Aiw) = 
y^. fi(aiw). As the fi are linearly independent over C, it follows that AiW = aiW for all i as claimed. □ 

3. Tensor Constructions on V and on Ad 

Before turning to general tensor constructions, we start with the (simpler) case of the action of G on 
symmetric powers of V to ease the exposition. The material in this section is mostly known to specialists 
but is here to set notations and clarify what follows. 
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3.1. Actions of G and Lie(G) on C[X\, . . . , X n ] 

Definition 11. Let V denote a C-vector space of dimension n and let G C GL(V) be a linear algebraic 
group. Let (Yi) be a basis ofV. Let a £ G be an automorphism ofV and M = (m^j) £ GL(n,C) its matrix 
in the basis (Yi). One defines an action (right translation) of a on C[X±, . . . ,X n ] by 

n 

cr(Xj) = J~]mj,j ■ Xi. 

i=l 

Given a homogeneous polynomial P £ C[X\, . . . , X n \— m of degree m, the action of a on P is defined by 

o-(P)(X u ...,X n ):= P(<r(*i), ■ ■ • . <r(X n )). 

We identify P with its (column) vector wp of coefficients in the symmetric power Sym m (V). Naturally, 
there exists a matrix, denoted by Sym m (M), such that 

v<r(p) = Sym m (M) ■ v P . 

This way, the constructor Sym m is a group morphism (from GL(V) to GL(Sym m (V)), called the m-th 
symmetric power in the sense of Lie groups. 

Definition 12. In the above notations, let again g = Lie(G). To a matrix N := (n.ij) £ q C Mat(n, C), we 
associate the derivation 

j=l \i=l / J 

The action of N G g on elements C[Xi, . . . ,X n ] is defined as the action of the derivation Dn on the 
polynomials P € C[X\, . . . , X n j. 

This may be understood using the e- formalism. Let e be a variable subject to e 2 = ; we know ([vdPS03], 
chap. 1, or [MS02], part 5) that N £ g if and only if Id + eN satisfies the group equations (modulo e 2 ). If 
we set <7 e := Id + eN, the above action (and the Taylor formula) induce cr £ (P) = P + eDn{P), hence this 
definition. As above, this defines a matrix sym m (iV) such that 

v Dn (P) =sym m (iV) -v P . 

We say that sym m (./V) is the m-th symmetric power of N in the sense of Lie algebras. 

Remark 13. The calculation following definition 12 shows that, for B G Mat(n, k), an alternative charac- 
terization ofsym m (B) is 

Sym m (Id + sB) = Idjv + £ sym m {B). 

Similarly, we consider a linear differential system [A] : Y' = A ■ Y with A £ Mat(n, k) and let U be a 
fundamental solution matrix; using again a variable e satisfying e 2 = and using section 3.3 below to define 
evaluation of U (if needed), we have U(x + e) = U(x) + eU'{x) = (Id + eA)U{x) so Sym m (?7)(a: + e) = 
(Id + e sym m (A))Sym m (C/)(x). It follows that 

Sym Tn (J7)' = sym m (A) ■ Sym m ([/). 

Lemma 14. Consider the differential system [A] : Y' = AY with A £ Mat(n, k) and the change of variable 
Y = PZ with P £ GL(n, k). Let B = P[A] := p- 1 ■ (A ■ P - P'). Then, we have 

S ym m (P[A}) = S ym m (B) = Sym m (P)[ S ym™(A)]. 

Proof. Let U and V be fundamental solution matrices of [A] and [B] respectively. As U = P ■ V and Sym m 
is a group morphism, we have Sym m (C/) = Sym m (P) • Sym m (V) hence the result. □ 
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3.2. Tensor Constructions and Differential Modules 

As above, V denotes the solution space of Y' = AY in a Picard-Vessiot extension K of k with the 
differential Galois group G acting on V. 

Definition 15. A tensor construction Const(y) on the G-module V is a vector space obtained from V by 
finite iteration of ' <S>, ©, * (dual), symmetric powers Sym m and exterior powers A r . 

Remark 16. In representation theory, one also considers subspaces and quotients. Here, these "direct" 
tensor constructions (GL„ -modules) will be enough for our purpose. The proofs of the main results of this 
paper rely on Chevalley's theorem ([Hum75] thm 11.2, [Bor91] thm 5.1, [Spr98] thm 5.5.3) which states that 
there exists a tensor construction Const(y) and a line D C Const(T^) such that G is exactly the stabilizer 
of D in GL(V) ; the fact that the tensor constructions used here are enough can be seen by inspecting the 
proofs of Chevalley 's theorem; for example in [Ber86] page 67, it is shown that there is a finite dimensional 
W C C[End(V)] ~ Sym(V (8 V*) such that D — A r (W) (and hence lies in a tensor construction) . 

The way G acts on tensor constructions can be found e.g. in [Hum75, FH91] or chapter 2 of [vdPS03] ; in 
particular, given a basis Y = (Yi)i = \... n of V and a a tensor construction Const(V^) on V, one may construct 
a canonical basis Const(Y) of Const(V). Given g € G with matrix M, the matrix Const(M) is defined (as in 
the above subsection) as the matrix of the action of g on Const(y) relatively to this basis. The constructor 
Const is then clearly a group morphism from G to GL(N,C) (with N = dim(Const(V r ))). The associated 
"Lie algebra" constructor const may be (as above) defined by the identity Const(Id+eA r ) = Id^ + e const(TV), 
i.e. by the way a derivation Dn (in the above notations of section 3.1) acts on a construction. This makes 
const a vector space morphism. 

Lemma 17. Let A — £\ /jAj be a Wei-Norman decomposition of A. Let Const(») be a tensor construction. 
Then const(A) = J^. / i const(A i ) is a Wei-Norman decomposition of const(j4). 

Proof. Follows from linear independence (over C) of the fa and the fact that const acts linearly. □ 

Let now A4 := (k n , Va = d — A) be the differential module. As shown in chapter 2 of [vdPS03], the 
matrix of the action of Va on Const(A^) is const(A). In particular, in terms of a fundamental solution 
matrix U G GL(n, K), we have 

Const(C/)' = cotlSt(A) • Const([/). 

3.3. A good fundamental solution matrix and its evaluation 

In the rest of this paper, we'll choose a convenient fundamental matrix U. If k = C(z) then we choose 
some ordinary point zq £ C of [A] and let U denote the fundamental matrix of [A] in C[[z — zo]] such that 
U(z ) =Id. 

For a more general differential field, we can do more or less the same (we follow here the treatment in 
[BroOl], section 3). Let k denote a finitely generated differential extension of Q. The n 2 coefficients of A 
then lie in a differential field k which is a finitely generated differential extension of Q. By Seidenberg's 
Embedding Theorem [Sei58, Sei69], any such field is isomorphic to a differential field T of meromorphic 
functions on an open region of C. 

So there exist infinitely many zq in C such that all coefficients of A can be seen to be analytic in an 
open neighborhood of zq (containing zq). Such a zq will be called an ordinary point of [A]. We may then 
apply Cauchy's theorem to construct a local fundamental solution matrix U £ C[[z — zq]] normalized by 
U{z ) - Id. 
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Once U is chosen, we first note that det(f7)(zo) = 1- Next we recall that, for a tensor construction 
Const, the entries of Const(J7) are polynomials (with coefficients in C) in the entries of U and 1/ det(J7). 
As &et{U)(z ) = 1, the inverse of det(U) lies in C[[z — z ]} so that the entries of Const(J7) arc polynomials 
in elements of C\[z — zq]]. As a consequence, the constructor Const commutes with evaluation at zq and, as 
Const is a group morphism from GL(n, C) to GL(N, C), 

Const(J7)(z ) = Const([>(zo)) = Const (Id„) = ld N 

In the rest of the paper, the reader may thus think of k as being C(x) (or C{{x})) and then use this 
embedding to obtain the results for general k. 

Remark 18. Lemma 17 shows that, when one performs tensor constructions (in the sense of definition 15) 
on a differential system [A] , no apparent singularity may appear ( this is not true if one performs constructions 
on operators). Indeed, the coefficients o/const(A) will be linear combinations of the coefficients of A. Hence, 
if we choose an ordinary point Zq of [A], it remains an ordinary point of any tensor constructions const(A). 

3.4- Invariants and Semi- Invariants of G and of Ai 

Definition 19. Let Const(V) denote a tensor construction on V. An element I £ Const (V) is called a 
semi-invariant of G if, for all g € G, g(I) = Xgl where x '■ G — > C* is a character of G. 
We say that I is an invariant of G if, for all g £ G, g(I) = I . 

For / in k, we write exp(J /) for a solution of y' = fy. A solution of Y 1 = AY is called rational if Y £ k n ; 
it is called exponential if there exist f £ k and F £ k n such that Y = cxp( J f) ■ F. 

Definition 20. A rational (resp. exponential) solution of some construction Y' = const(A)Y will be called 
an invariant (resp. semi-invariant 3 ,) of the differential module M.. 

As in section 3.3, pick an ordinary point zq £ C and let U denote the fundamental matrix of [A] in 
C[[z — zq]] such that U(zq) = Id. To an clement of V, represented by its vector v of constants on the 
basis given by the fundamental matrix U, we associate the element <j> v := Uv of A4 <g> K ; this gives the 
isomorphism V (£> K — > M. ® K. For a tensor construction Const on V, to a v £ Const(V r ) we associate the 
element <f> v := Const(C/).ii £ Const(A / t) <£> K. For self-containedness, we recall the following standard result. 

Lemma 21. In the above notations, I is an invariant in Const(T^) if and only if <f>j is a rational solution 
of [const(-A)], i.e. <j>i = t (Xi, ■ ■ ■ , A m ) with A, £ k, for all i. 

Similarly, I is a semi-invariant in Const(V) if and only if (pi is an exponential solution of [const(A)]), i.e. 
4>i = cxp(J /) '(Ai, • • ■ , A m ) with f £ k and A, £ k, for all i. 

Proof. Let g £ G. It acts on U by g(U) = U.M g . Given an element I £ V represented by its vector v of 
coefficients in C, the coefficient vector of g(I) is M g .v. If J is an invariant in Const(y), then Const(M s ).w = v 
so g(4>i) = 4>i, hence f>i £ k n . Similarly, if / is a semi-invariant, then Const(il/ 5 ).w = Xg v an d g{4>i) — Xg4>i- 
Pick a non-zero coordinate, say /i, of <j)j. As g(fi) — Xgfi f° r a U g £ G, we have £ k. Now, for any 

other coordinate ft of <fii, g(fi/ fi) = fij fi so f%lf\ £ k. □ 

This lemma gives the standard dictionary between rational (resp. exponential) solutions of Y 1 = 
const(A)y and invariants (resp. semi-invariants) in Const(V r ). 



3 If (j>i = exp(f f) f (Ai, • • • , A m ) with A; S k, then it is the vector '(Ai, ■ ■ ■ , A m ) which lies in the module and could be called 
a semi-invariant of the module, so our notation is a bit abusive but coherent with the dictionary below. 
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Dictionnary 22. If v G C N is the vector of coefficients of an invariant I G Const(V^), the corresponding 
invariant of the module (rational solution ofY' = const(A)Y ) is </>/ = Gonst{U).v. 

Conversely, given an invariant <fi of the module (rational solution ofY' = consi(A)Y), we use the evaluation 
technique from section 3.3 : as U{ztf) = Id 7 we will have <f> = Const(J7).w where v = 4>(zo) is the coefficient 
vector of the invariant I G Const(V) associated to (f> (because Const(Id) is again the identity). 

A similar dictionary is used e.g. in [Bro01],[CS99], [vHW97]. We may now introduce a key definition for 
our reduction criteria. 

Definition 23. Let [A] : Y 1 = AY with A G Mat(n, k) be a linear differential system. Let Ai := (k n , V = 
d — A) be the associated differential module with basis {e^}. 

An invariant (f> G Const(yVJ) v of the module (a rational solution (f> = *(Ai, . . . , Ajy) of <ft' = const(A)(f>) 
will be said to have constant coefficients when A,; G C for all i = 1 , . . . , N . 

A semi-invariant <p G Const (.M)<8>i^ of the module (an exponential solution </>/ = exp(J /)• (Ai, • • • , Ajv), 
with f G k and Xi G k, of <f> = consi(A)(f>) will be said to have constant coefficients if f can be chosen so 
that Xi G C for all i = 1, . . . , N. 

4. Characterization of Reduced Forms via their Semi-Invariants 

In this section, we give our main result, a simple criterion to constructively characterize a system in 
reduced form. 

Theorem 1. Let k denote a finitely generated differential extension of Q. Let [A] : Y' = AY with A G 
Mat(n, k) be a linear differential system. Let Ai := (k n , V = d — A) be the associated differential module. 
The system [A] is in reduced form if and only if any semi-invariant of the module Ai has constant coefficients 
(in the sense of definitions 20 and 23). 

Proof. As in section 3.3, we fix a fundamental solution matrix U of Y' = AY at an ordinary point zq such 
that U(zq) — Id„. This choice of basis for V induces a matrix representation of G and q in GL{V). 

Assume that Y 1 = AY is in reduced form. We therefore have a Wei-Norman decomposition A = ^ fiAi 
(fi in fc, linearly independent over C) where each Ai G g. 

Let J be a semi-invariant in a tensor construction Const(V). We let v :=* (a\,--- ,ajv) G C N be the 
coordinates of / in Const(V) relatively to the basis induced by U. We let, as in dictionary 22, <pj := 
Const(J7) • v. As / is a semi- invariant, the associated exponential solution of Y' = consi(A)Y is <pj = 
exp(J /) •* (Ai, • • • , Xn) with / and the A.; in k. Of course, as Const([/)(2;o) = Idw, we have 0/(zo) = v. 

We will show that the Ai may in fact be forced to be constants. Recall that const(A) decomposes as 
cortst(A) = fiConst(Ai). As Ai G q and v is the coordinate vector of the semi- invariant /, there exist 
constants Ci such that const(Ai)v = CiV. Define tp := exp (J Cifij and choose the integration constant 
such that <p(zq) = I. Let F be the vector in K" defined by F := tp ■ v. As the ai are constant, we have 

f' = <p r -v = (p (y^ c */*) v - 

Now 

const(A) ■ F = tp I ^2 fiConst(Ai)v = tp f ^ cjA v. 
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It follows that F' = const(A) ■ F. As <p(zo) = 1, we have F(zq) = v. So F and <pi are two (exponential) 
solutions of Y' = const(A)Y with the same initial condition 4>i(zo) = F(zq) = v at an ordinary point Zq, 
hence (j>i = F and 4>i indeed has constant coefficients. 

We now prove the converse. Assume that any semi-invariant of the system Y' = AY has constant 
coefficients (in the sense of definition 23). By Chcvalley's theorem (see remark 16), there exist a tensor 
construction Const(T^) and an clement / G Const(l/) such that the algebraic group G is exactly the set of 
automorphisms of V which leave the line generated by / stable. If we let v denote the coordinates of / in 
Const(V^), then 

q = {M G Mat(n, C) \3c G C, const(M )v = cv} . 

As above, let <\>i :— Const(J7) ■ v be the associated exponential solution of Y 1 = const(A)Y . Note that 
4>i( z o) = v. By hypothesis, there exist an exponential ip (i.e. ip G K and 3/ € k such that tp'/ip = f) and a 
vector w <G C N of constants such that <pj = tpw. If we choose the integration constant such that ip(z ) = 1, 
we have w = v (evaluate in zq). 

Let A = J2i fiM be a Wei- Norman decomposition of A. Then const(A) = J2i fi cons K^i) i s a Wei-Norman 
decomposition of const(A) (lemma 17). Now (f>i is a solution of Y' = const(A)Y ; we compute 

= 0j - const(A)0 7 = ip (^fv - /iC0tist(Ai)u"j . 

So v is a constant eigenvector of fiCOXisi^Ai). By lemmas 17 and 10, it follows that v is an eigenvector of 
each const(Ai) which in turn implies that, for all i, we have Ai € Q. Hence A is in reduced form. □ 

A vector v £ Const(V) is called a Chevalley semi-invariant for the group G when 

G = {M G GL(n, C) \3c G C, Const(M)u = cv} . 
The above proof shows the following finer corollary : 

Corollary 24. In the above notations, the system [A] is in reduced form if and only if there is a Chevalley 
semi-invariant v in a tensor construction Const(y) of G such that (f> v := Const(C/)w has constant coefficients. 

Another consequence of the proof of this theorem is that it may be adapted to a result for partial 
reduction. 

Corollary 25. Consider a differential system Y' = AY with A G Mat(n, k). Let G be the differential Galois 
group and H D G be some group containing G with Lie algebra f) = Lie(H). Then A G f)(fc) if and only if, 
for any semi-invariant I of G which is also a semi-invariant of H , its image <pi has constant coefficients. 

Proof. The implication is proved like in the theorem. For the converse, replace / in the above proof by a 
Chevalley semi-invariant defining H (instead of G). □ 

We may illustrate this result with a very simple example. If G C H := SL(n,C), then the determinant 
det of a fundamental solution matrix is an invariant - and this is the only invariant of G which is an invariant 
of H . The equation that it satisfies is det' = Tr{A) det. Here, f) is the set of matrices with zero trace; we 
see that A is in t)(k) if and only if det' = 0, i.e. if det is constant. 
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5. A Reduction Procedure when G is Reductive and Unimodular 

When G is reductive and unimodular, we may refine our reduction criterion by reducing it to invariants 
in Sym(nV). Furthermore, we give a reduction procedure in this case (which also gives a constructive proof 
of the Kolchin-Kovacic reduction theorem). In what follows, we assume that G is reductive and unimodular; 
the term "Invariant" will thus refer specifically to an invariant in Sym(ny) (and an invariant of the module 
will be an invariant in Sym(n.M) v ). 

5.1. A Characterization of Reduced Forms via Invariants when G is Reductive and Unimodular 

When G is reductive, we start by recalling a useful variant of Chevalley's theorem (see also [Bor91], 
chap. II §5.5, page 92) 

Lemma 26. Let G C GL(n, C) be a reductive linear algebraic group acting on an n-dimensional vector space 
V . There exists an invariant I in a tensor construction on V such that G = {g E GL(n, C) \ g(I) = I}. 

Proof. By Chevalley's theorem (see remark 16), there exists a tensor construction W on V and a one 
dimension subspace V = C.v of W such that G — {g E GL(n,C)\g(v) = Xg v iXg € C*}. As G acts 
completely reducibly, V admits a G-stable supplement V in W so that the dual W* decomposes as a direct 
sum W* — V* © V* of G-spaces and G acts on v* as g(v*) = 1/xgV* (see e.g. [Wei95] lemmas 15 and 16 
and their proofs). Then, / := v <g> v* is an invariant of G in W <£> W* and it is easily seen that G = {g E 
GL(n,C) \g (/)=/}. □ 

Again, such an invariant / E Const(y) giving G = {g E GL(n, C) \ g(I) = 1} will be called a Chevalley 
invariant for G. For a reductive group G, the proof of theorem 1 can be straightforwardly adapted (and we 
let the readers adapt its two corollaries to the case of invariants too) . 

Proposition 27. Let k denote a finitely generated differential extension of Q. Let [A] : Y' = AY with 
A E Mat(n, k) be a linear differential system. Let M. := (fc" , V = d — A) be the associated differential 
module. We assume that the differential Galois group Gal([A\) of [A] is reductive. 

The system [A] is in reduced form if and only if any invariant of the module M. has constant coefficients ( in 
the sense of definitions 20 and 23) 

Proof. The implication is proved as in theorem 1 (invariants are semi-invariants) as follows. Let I be an 
invariant in a tensor construction Const(y). We let v :=' (ax,-- - ,av) E C N be the coordinates of / 
in Const(V r ) relatively to the basis induced by U (from section 3.3). We let, as in dictionary 22, 0j := 
Const(J7) - v. As I is an invariant, the associated rational solution of Y' = const(A)Y is <f>i =' (Ai, • • ■ , Ajv) 
with the Ai in k. Of course, as Const([/)(zo) = Idv, we have <j>i(zo) = v. Recall that const(A) decomposes 
as const(j4) = fiCOX\st(Ai). As A; E q and v is the coordinate vector of the invariant /, we have 
const(Ai)v = so const(A)u = £) 4 /jCOtist(Ai)v = 0. As v' = 0, we see that v' = const(A)v. Now, as 
4>i( z o) = v i we see t na t v and 4>i are two (rational) solutions of Y' = consi(A)Y with the same initial 
condition at an ordinary point Zq, hence <pi = v which has constant coefficients. 

To see the converse, just use the above lemma 26 to replace " semi- invariant" by " invariant" and " exponential" 
by "rational" in the proof of the converse in theorem 1. □ 

We now assume in the sequel that G is reductive and unimodular and prove a finer version which restricts 
the type of constructions used then. For this we recall some facts. 

Let Sym(n^) denote the ring of tensors of Sym(nV) := Sym(V^ © • • • © V) (n copies of V) which are 
invariant under the action of G. As G is reductive this ring is finitely generated. Let Ji, • • • , I p be a system 
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of generators of this ring. As G is reductive and unimodular we can use [Com96] (or [CS99]) and convert this 
basis of Sym(nF) G to a basis Pi,- - ■ ,P p of the ideal J of algebraic relations with coefficients in k satisfied 
by the elements of a fundamental matrix of solutions U of the differential system. We can define an action 
of G on the polynomials of by setting aP(Xij) = P(XijM a ) where M a is defined by a(U) = UM a as 

in section 3.1. 

With this action, the polynomials Pj are invariant (because the tensors Ij are). The Galois group G is 
then defined as the set of automorphisms of V that preserve the ideal J ([CS99] or [vdPS03]) and this is 
equivalent to say that G is defined by 

G = {a e GL(V) | a(I k ) = h,l<k< p}. 

We also deduce (as in section 3.1) that the Lie algebra Lie(G) of G is defined by 

Lie(G) = {g e Ql(V) \ D g (I k ) = 0, 1 < k < p}. 

Theorem 2. Let k denote a finitely generated differential extension of Q. Let [A] : Y' = AY with A G 
Mat(n, k) be a linear differential system. Let M. := {k n , V = d — A) be the associated differential module. 
We assume that the differential Galois group Gal([A]) of [A] is reductive and unimodular. 
The system [A] is in reduced form if and only if any invariant <j> € Sym(n,.M) v of the module has constant 
coefficients (in the sense of definitions 20 and 23). 

Proof. Assume that [A] is in reduced form. By proposition 27, all invariants have constant coefficients so we 
only need to prove the converse implication. Assume that, for any construction Const (M.) = Sym m (n.M), 
any invariant </> € Const(.M) v has constant coefficients (in the sense of definition 23). Using section 3.3, we 
choose an ordinary point Zq and a fundamental solution matrix U verifying U(zq) = Id. The notation being 
as in the previous paragraph, we take the columns of U as a basis of V and express the coordinates Vk of 
the generating invariants Ik relatively to the induced basis in constfe(V r ) := sym mk (nV). We thus have 

= {N € Mat(n,C) | const fc (iV) • v k = 0, 1 < k < p}. 

Let (f>k be the image of Ik in Constfc(?i.M) v given by dictionary 22, i.e 4>k = Constfc(J7) • Vk- As (j>k is constant 
by assumption (and Constfc([/)(2;o) is the identity), the evaluation of this relation at zq gives <f>k = Vk- 

Now, as <f)' k = constk(A) ■ <f)k and constfc(A) = J2l=i fa ' constfc(Ai) is a Wei-Norman decomposition of 
cortstfc(A) (lemma 17), we have Yll=i fa corls tfc(^L) ■ Vk = 0. And, as the fi are linearly independent over C, 
it follows that constfc(Ai) • Vk = for alH = 1, . . . , r and for all k = 1, . . . ,p which proves that Ai G g for 
i = 1 , . . . , r and the system is in reduced form. □ 

Remark 28. This proof sheds light on another useful feature (known to differential geometers, see e.g. 
[DNF85], Chap. 6 §25 lemma p. 238). If the system is in reduced form, the above proof shows that the 
canonical local solution matrix U satisfies Vk = Constk(U)vk for k = 1, . . . ,p; as G = {g € GL(V) \ g{Ik) = 
Ik, 1 < k < p} (see the discussion preceding the theorem), this means that the canonical local solution matrix 
U satisfies the equations of the group : U £ G(K). 

There remains a minor point in order to make this effective. Assume that we start from a matrix 
A G Mat(n, k). If we compute a reduction matrix P, it will have coefficients in some algebraic extension ko 
of k. So, to apply our theorem, we would need to compute solutions of the Sym m (nAi) <g) ko in kQ and we 
would like to avoid that. This is provided by this simple but useful lemma. 
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Lemma 29. Let [A] : Y' = AY denote a linear differential system with A £ Mat(n, k). Assume that we 
know an algebraic extension ko of k and a matrix P £ GL(n, ko) such that, for any invariant <fi £ Const(.M) v 
(i.e. <fi is a solution in k N of Const (A4)), Const(P) -1 ^ has constant coefficients. Then, for any (algebraic) 
solution (f> of any (Const(.M) ® k), Const(P) _1 also has constant coefficients. 
In other words, P[A] is in reduced form and P is a reduction matrix. 

Proof. Take an algebraic solution <j> of some construction Const (i.e. <fi £ (Const(A^) <8> fc) v ) and choose 
an entry / of this solution. The coefficients of the minimal polynomial of / are given as entries in k of a 
solution of another tensor construction : as shown e.g. by Singer and Ulmer [SU93, SU97], coefficients of 
minimal polynomials of algebraic solutions of linear differential equations are found as invariants in symmetric 
powers. As P _1 maps the latter to constants, the image of / under P~ J will have a minimal polynomial 
whose coefficients are constants. As C is algebraically closed, this in turn implies that the image of / is in 
C, hence the result. □ 

5.2. The Compoint-Singer Procedure for Computing G when G is Reductive 

We now assume that k is an algebraic extension of C(x). Theorem 1.1 of Compoint in [Com96, Com98] 
shows that if a Picard-Vessiot extension K has a reductive unimodular Galois group G, then the generators of 
the ideal of relations among solutions (and their derivatives) are obtained from the generators of Sym(nA^) v . 
From van Hoeij and Weil [vHW97] we know that, for a given m, the generators of degree m of Sym(nA4) v 
may be constructed directly from system [A]. So computing G is reduced to computing a bound on the 
degree of the generators of the ring of invariants of G. The work of Compoint and Singer [CS99] shows how 
to compute such a bound and culminates in this result: 

Proposition 30 (Compoint, Singer [CS99]). Let k denote an algebraic extension of C(x), where C is an 
algebraically closed computable field. Let [A] : Y' = AY be a linear differential system with A £ Mat(n, k) and 
a reductive unimodular Galois group. Then one can compute in a finite number of steps a basis (f>i, . . . , <j> r £ 
k[Y itj ] /Sym(nX) v . 

From an implementation viewpoint, this procedure is far from being efficient, but it gives a (yet somewhat 
theoretical) starting point for our reduction algorithm below. 

5.3. A Reduction Procedure for Reductive Unimodular Differential Systems 

We assume, as above, that k C C(z) and that G := Gal ([A]) is reductive and unimodular. Continuing 
with our notations, we consider U £ GL(n , C((z — Zo))) a fundamental solution matrix of [A] satisfying the 
initial condition U(zq) = Id at an ordinary point zq of [A]. At this point, we have a simple procedure for 
computing a reduced form. 

Reduction procedure 

1. Choose a regular point zq in C and a canonical local solution matrix U £ GL{n , C((z — zq))) satisfying 
U(z ) = Id. 

2. Using the Compoint-Singer procedure, compute a set of invariants <pi, . . . , <fi r in Sym(nA^) v and the 
set of invariants I { := 4>i(z Q ) of G in Sym m *(W) such that G = {g £ GL(V)\Vi,g(Ii) = /,}. 

3. Pick a matrix P with indeterminate coefficients. Define a system (S) of polynomial equations 

(S) : Vi = 1 • • • , r, Sym m< (P ® ■ ■ ■ ® P) ■ J 4 - ^ = 0, det(P) ^ 0. 
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4. Using any algorithm for polynomial system solving (Groebner bases, triangular sets), find a solution 
in GL(n, k) of (S). The resulting P[A] is in reduced form. 

Theorem 3. The above algorithm is correct and computes a reduced form for any system with a reductive 
unimodular Galois group. 

Proof. Recall that an algorithm is correct when each step is well defined, it terminates in finite time, and 
the result is what was expected. 

The fact that Step 1 makes sense follows from section 3.3. That Step 2 is well defined follows from [CS99] 
(using [Com96, Com98] and [vHW97]). Step 3 is well defined if we can guarantee that system (5) is always 
consistent. But, by construction, the local fundamental solution matrix U satisfies (5). As (S) has this 

2 — 2 

solution in K n , by Hilbcrt's Nullstcllensatz, it has a solution in k n (which can be obtained by Groebner 
bases [CLO07] or triangular sets (e.g. [Hub03])). Now our theorem 2 and lemma 29 show that the resulting 
P[A) is in reduced form. Indeed, P sends the & to invariants Z; with constant coefficients and, by assumption, 
any invariant is polynomial in the 4>i and is hence mapped to one with invariant coefficients. □ 

Remark 31. This algorithm computes a reduction matrix with coefficients in k. However, when G is 
connected, the result of Kolchin and Kovacic (proposition 3) shows that there should exist a solution P with 
coefficients in k. This is where the C\ condition on k enters. If the equations defining P are quadrics over 
k, then there exists an algorithm to find a solution in k ([CvH06]). In general, though the proof of Tsen's 
theorem can be made constructive, we do not know of an implemented algorithm to determine a rational point 
on an algebraic manifold on k in full generality. So the above algorithm, though correct and complete, does 
not fulfill the Kovacic program completely in the sense that it does not yet address the rationality problem 
"find P € GL(n,k)" in the case when G is connected. 

Note that, apart from the rationality issue, this theorem gives us a constructive proof for the classical 
Kolchin-Kovacic reduction theorem (proposition 3 above) for reductive unimodular groups. We know (e.g. 
from section 5 of [MS02]) that, if A <G fy(k), then there exists a fundamental solution matrix in H(K) so that 
(S) will have a solution in H(K) as in the above proof. 

To obtain a similar proof of the reduction theorem for general groups (dropping the C\ assumption on 
the fields and the rationality issue), pick a Chevalley semi-invariant (f>i exponential solution of some const(A) 
as in the proof of theorem 1 and write (j>i = exp(J ip)V with exp(J (p)(zo) = 1 and V G k N ; letting 
/ := V(zq), the same reasoning shows that, for a P with unknown coefficients, the system V = Const(P)./ 
has a solution P in GL(k) (because U satisfies these equations so they are consistent) and the proof of 
theorem 1 then shows that P[A] is in reduced form. If the group is connected and we have a reduced form 
with coefficients in k then to find one with coefficients in k is a descent question (see e.g. [CvdPWIO, CW04]). 

We now turn to the structure of the solutions of the polynomial system (S) of the reduction procedure. 
We first note a simple lemma. 

Lemma 32. Let k be as above and A £ Mat(n,k). Let kg be an algebraic extension of k and B £ Mat(n, ko) 
such that the system Y' = BY is a reduced form of [A] . 

1. The differential Galois group of [B] is connected. 

2. // there is a reduction matrix P £ GL(n, k) such that B = P[A], then ko = k and Gal([A]) is connected. 

Proof. We assume that kg is the coefficient field of B, i.e. the smallest field containing all coefficients of 
B. Let K be a Picard-Vessiot extension of fco associated to [B] and G = Gal(K/ko) = Gal([B\) be the 
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differential Galois group with Lie algebra g. As B £ fl(fco), Kolchin's theorem (proposition 1.31 in [vdPS03], 
take H = G°) shows that G C H := G°, hence G = G° and G is connected. If now B = P[A] with 
P £ GL(n, k), then the coefficient field of B is inside k and Gal([B\) = Gal([A]) so (2) follows from (1). □ 

Note that if A; is a C\ field then this lemma and the Kovacic-Kolchin reduction (proposition 3) show that 
Gal([A\) is connected if and only if there exists a reduction matrix P in GL(n, k). 

Proposition 33. Consider the polynomial system (S) of the reduction procedure. Let ko := K G denote the 
fixed field of the connected component G° of G in K. 

1. If k is a C\ field, then there exists a solution to (S) in GL(n,ko). Furthermore, ko is the smallest field 
over which a solution can be found. 

2. Suppose we have a solution P £ GL(n, k) of (S). A matrix P ■ M is another solution of (S) if and only 
ifM£G(k). 

3. A matrix P £ GL(n, k) is a solution of (S) only if P' ■ P^ 1 = A + N with N £ g(k). 

Proof. View the differential system as having coefficients in ko- Then its Galois group is G° and is connected. 
So, by proposition 3 of Kolchin and Kovacic, there exists a reduction matrix in GL(n, ko), hence Part 1. The 
fact that fco is minimal follows from the previous lemma. 

For part 2, we have Const(PM)./ = <f> = Const(P)./ so, as Const is a group morphism, Const(M)./ = / 
for all invariants defining G. This proves that M £ G(k). Part 3 is the derivative of this relation. Let 
B := P' ■ P- 1 We have <j>' = const(A)^ and 

4>' = Const(P)7 = const(B)Const(P)7 = const(B)(j) 

so, as the constructor const is linear, const(P — A)<fi = and B is given as the solution of a linear (non- 
differential) system. □ 

Note that even the reduction matrix P has coefficients in a big algebraic extension, the reduced form 
P[A] may still be defined over a small field. For example, consider a system whose solutions are all algebraic. 
Then P is just a fundamental solution matrix (algebraic) and P[A] is the zero matrix (which has coefficients 
in a rather smaller field). 



6. Examples 

6.1. A Dihedral Case in GL(2) 




A variant of Kovacic's algorithm shows that the group is a central extension of the infinite dihedral 
group. Indeed, it has an invariant in Sym 2 (A 2 (A / l)) (with value x) and and invariant in Sym 2 (A'l) with 
value (— 1,0,2;). If we choose the evaluation point Zq = 1, the corresponding invariants are I\ = 1 and 
I2 = (—1, 0, 1). Following our procedure, the reduction matrix P should satisfy 

jz - (pi,lP2,2 -Pl,2P2,lf ,2p M p 2 ,l - 2pi, 2 P2,2, -1 + Pl,l 2 - Pl,2 2 , X + p 2 ,l 2 ~ P2.2} 
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A triangular set computation gives simpler equations (pi^l + p\^ 2 ,x + £>2,i 2 , 2*2,2) hence the reduction 
matrix 

i 



and 



P := 



P[A] 



i\fx 

1 

1 



If we want to recover the more standard diagonal representation, we conjugate again by 

1 -1 



Po 



1 1 



which turns our invariant [—1, 0, 1] into (0, 1, 0) and 

(P.P 2 )[A] = v 




which is obviously a reduced form. Note that, in this example, it is not possible to obtain a reduced form 
without extending the ground field. 

6.2. Example: the Reduction Method when Lie([A]) = so(3) 
The special orthogonal group SO(3) is defined as 

SO(3) := { g £ GL 3 (C) : g € stab(Y 1 2 + F 2 2 + Y 3 2 ) et det(g) = 1}. 

We consider A <S Mat(3, k) and show how to put [A] : Y' = AY in reduced form with our procedure when its 
Galois group is SO(3). Note that as SO(3) is conjugate to a symmetric square of SL(2), this is an instance 
of the old problem of solving linear differential equations in terms of equations of lower order (see [Sin85] , 
[Per02][vH07, vHvdP06], [Ngu09, NvdPIO]). 

Testing whether G C S0 3 (C) or not is achieved in three steps (see [SU93, vHRUW99]): 

1. To check whether G C SL(3 , C) :={je GL(C, 3) : det(g) = 1}, check whether there exists w £ k 
such that w' = Tv(A)w. If yes, apply to A the gauge transformation: 



P 



uj 
1 
1 



so that Tr(P[A]) = 0, i.e. P € st(3). Remark that if U is a fundamental solution matrix of [A], then 
V := P~ 1 U will be a solution matrix of Z' = P[A]Z such that det(V) is constant, in the spirit of our 
criterion. 

2. Using any factorization algorithm, check whether [A] is irreducible (if it is reducible, use [AMW12]). 

3. Find an invariant of [A] of degree 2, i.e. test whether [sym 2 (A)] has rational solutions. If it does not, 
then q <f_ so(3). Otherwise, we have a solution Y — t (fx, f 2 , fs, fi, fs, fe) S k 6 of Y 1 = sym 2 (A)Y . 
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Now, to achieve reduction, we find a gauge transformation Q which transforms Y to the form *(1, 0, 0, 1, 0, 1). 
Indeed, if we let U be a solution matrix satisfying U(zq) = Id„ then Y = Sym 2 (J7) •/ and, by Gauss reduction 
of quadratic forms, I is conjugate to (1,0, 0, 1, 0, 1) i.e. X\ + X| + X%. 
We apply this to the system given by the matrix 



.4 



2x 2 -2x+l 
x(-l+x 2 ) 
x(2x-l) 
' {x + l){-l+x 2 
x 2 (x-l) 
x+1 



5x-3x 3 +2x 5 -1+x 2 -x 4 - 
( x — 1 ) X 
x 5 -x 4 -x 3 +x 2 +4x-l 



-1+x 2 

-x(x — 1)(1 — x 2 



x A ) 



2x^ -3x" +x+2 
x 2 (x — l) 2 
x 4 -2x 3 +2x 2 + l 
x(x+l)(x-l) 2 
5 -2x 4 +x 3 +2x-l 
x(x — l) 



We check that [A] is irreducible (e.g. with algorithms in Maple). The system admits an invariant of 
degree 2, i.e. a rational solution of [sym 2 (^4)], namely I2 ■= — 



_ (3-x 2 +x 4 -2x) Y 2 
x 2 (x-iy 2 ^1 



(x 2 



J x(x+l)(x 



2x+2) v v 
JY7A1A2 



(2x - 2)X 1 X 3 



(x 2 -2x+2) 
(x + l) 2 (x-l) 



2 *!s=p. x x 

x+1 z 



x 2 (x ~ 1) 2 X 2 . (which can be found using [vHW97] or the 



Barkatou algorithm applied to sym (A)). Using standard Gauss reduction of quadratic forms, the invariant 



I2 becomes: 



(-3 + . 



2x) Z\ 



(x 2 - 2x + 2) 
3 — x 2 + x A — 2x 



■Z 2 



x 2 - 2x + 2 



-Z 2 . 



We now apply an algorithm for polynomial solutions of quadric equations (e.g. [CvH06]) which yields the 
polynomial change of variables Z = P2 ■ Y defined by 



P 2 := 



{x + l){x-l) 

-1 

1 



x 



x 2 - 2x + 2 
(x + l)(x-l) 
(x + l)(x-l) 2 





1 -x 
1 



So, finally, we combine this into transformation (Xl, X2, X3) = (Yi, Y2, Y^).P with 

x(x — 1) 



p 



1 

. 2 -l -i 

x(x-l) 



and we obtain P[A] 






X 


1 


— X 





X 2 


-1 


-x 2 






eso(3)(fc). 
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